Short-time dynamics and high-frequency rheology of core-shell suspensions with thin-shells 
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Short-time dynamics and high-frequency rheology for suspensions of non-overlapping core-shell particles 
with thin shells were analysed. In the thin-shell limit, the single-particle scattering coefficients were derived and 
shown to define a unique effective radius. This result was used to justify theoretically (in the thin-shell limit) 
the accuracy of the annulus approximation with the inner radius equal to the effective hydrodynamic radius 
of the core-shell particle. The two-particle virial expansion of the translational & rotational self-diffusion, 
sedimentation and viscosity was performed. The virial coefficients were evaluated and shown to be accurately 
approximated by the effective annulus model, in contrast to the imprecise effective hard sphere model. 
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I. INTRODUCTION 

Recently, there has been a growing interest in micro and 
nanogels, and other permeable particles, which can be used to 
carry drugs or proteins For such systems, density of 

polymer segments inside the core region is frequently much 
higher than in the outer part, and therefore they are often ap- 
proximated as core-shell particles. 

A core-shell particle consists of a solid core of radius a, 
and a surrounding permeable shell, with the inner and outer 
radii, a and b, respectively (see Fig. [TJ. The porous medium 
inside the shell is characterized by the uniform hydrodynamic 
penetration depth k _1 . It is assumed that the particles do not 
overlap, i.e. their centers cannot come closer to each other 
than 2b. 




FIG. 1: The core-shell particle: the solid core of radius a is sur- 
rounded by the uniformly permeable shell of radii a, b and the hy- 
drodynamic penetration depth 

In a recent paper [7], the precise multipole method was ap- 
plied to extensively analyze dynamics and rheology of core 
shell particles. High-frequency viscosity, short-time transla- 
tional and rotational diffusion, and sedimentation coefficient 
have been evaluated numerically for a very wide range of per- 
meabilities of both wide and thin shells, with volume fractions 
up to </> = 0.45. It was concluded that for wide shells with a 
low permeability (n(b — a) > 5), the core is practically invis- 
ible and particles behave as uniformly permeable. 
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In the oposite limit of thin shells, it was shown that the core- 
shell particles cannot be well-approximated by hard spheres 
of radius b, even for a tiny, hardly permeable shell. However, 
the dynamics of concentrated core-shell systems is well repro- 
duced, if a thin permeable shell is replaced by the shell of the 
pure fluid (the so-called annulus model [25]). The difference 
between the hard spheres and annulus particles of the same 
radii b is very large; for example, 20-35% for the diffusion and 
sedimentation, and 60% for the viscosity, if (b — a) /a = 0.05, 
ku = 10, and (f> — 0.45. The shell influence cannot be ne- 
glected, and the essential reason is that it prevents particles 
from overlapping. 

This result indicates that a more detailed analysis of 
the thin-shell limit is needed. For dilute systems of non- 
interacting core-shell particles, the single-particle transla- 
tional & rotational diffusion coefficients and intrinsic viscos- 
ity determine the corresponding hydrodynamic radii a e ff 18- 
[Till , which, in general, are only tenuously connected to the 
geometrical radius a of the core-shell particle [8|]. A thin- 
shell analysis of such systems has been recently performed in 
Ref . H] . The leading terms of a e ff have been given and shown 
to be the same for the intrinsic viscosity and the translational 
diffusion. 

In this paper, we extend the idea of Ref. |8f], and we in- 
vestigate the thin-shell limit for concentrated systems of the 
core-shell particles. We use the general concept of the hydro- 
dynamic radius model (HRM) [12], where the effective radius 
approximation is restricted only to the hydrodynamic interac- 
tion, with no modification of direct interactions. In this work, 
the direct interactions between core-shell particles are speci- 
fied by the no-overlap radius b. In this case, the HRM model is 
equivalent to the effective annulus model with the no-overlap 
radius of the outer shell b and the radius of the inner shell 
equal to the hydrodynamic radius a e g of the core-shell parti- 
cle. We will show that the effective annulus approximation 
is much more precise than the standard model of a hard-core 
particle with the effective radius a e ff . 

The paper is organized as follows. Sec.UJJcontains theoret- 
ical description of the core-shell suspensions. Sec. [TIT] is fo- 
cused on the hydrodynamics of the single core-shell particles 
with thin shells. In Sec. IIVI direct (no-overlap) interactions 
between core-shell particles are described within the annulus 
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model. In Sec. [V] the thin-shell limit of the two-particle virial 
coefficients is explicitly carried out for the translational & ro- 
tational self-diffusion, sedimentation and effective viscosity, 
and accuracy of the effective annulus model is examined. The 
conclusions are presented in Sec. [VI| 



II. THE CORE-SHELL SUSPENSIONS 

We consider a fluid of shear viscosity 770, with randomly 
distributed identical core-shell particles [9] defined in Sec. Q] 
and shown in Fig.Q] 

Geometry of the particle is characterized by the relative 
thickness of the porous shell, 



(1) 



relation can be written as a scattering series, 



The effect of the particle permeability is described by the ra- 
tio x of the particle radius a to the hydrodynamic screening 
length of the uniformly porous material inside the parti- 
cle, i.e. 



x = na. 



(2) 



All the distances are normalized by the core radius a. 

The fluid flow is characterized by the low-Reynolds- 
number, Re<<l. Outside the particles, the fluid velocity v 
and pressure p satisfy the Stokes equations fill, [Till . 



rj V 2 v(r) - Vp(r) = 
V-v(r) = 0, 



(3) 



and inside the particle shells, by the Brinkman-Debye-Biiche 
(BDB) equations id [Hi , 

?7o V 2 t>(r) — 770 k 2 [v(r) ~ Mj(r)] — Vp(r) = 

V-v(r) = 0. (4) 

The core and the skeleton of the particle i, centered at rj, move 
rigidly with the local velocity Uj(r) = Uj + fli x (r — rj, 
determined by their translational and rotational velocities U* 
and flj, respectively. The fluid velocity and the stress are con- 
tinuous across the particle surface, i.e. at |r — r^| = b. The 
fluid sticks to the core surface; that is, v = at |r — = a. 

The coupled Stokes-BDB prob lem defined above is solved 
by the multipole expansion lll7l [Hil . Many -particle mobil- 
ity coefficients are evaluated and averaged over the statistical 
equilibrium ensemble of configurations, resulting in the short- 
time self-diffusion, sedimentation and high-frequency viscos- 
ity coefficients, see e.g. Refs. |PtL fl9T[22ll for the general theo- 
retical scheme. The important issue discussed in this work is 
that the N-particle mobility coefficients /z i -(l...iV) are deter- 
mined in terms of the Green operators G(mn) (which depend 
on the distance between particles m and n), and the single- 
particle friction operators, Zo(k) and Zo(fc) (which depend 
on the internal structure and size of a particle k 112311 ). This 



fi, i Al...N)=6 ij fi (i) + fi (i) Z {i 



N 

E ft* 

k=l 



-(l-S ik )G(ik)Z (k)- 



(l-S kj )G(kj)Z (j)n (j). 

(5) 



or in short as 



1 



GZr 



G Z fi , 



(6) 



with fx denoting the single-particle mobility operator (deter- 
mined by the inverse of Zq). The single-particle friction oper- 
ators, Zq and Zo, determine the hydrodynamic force density 
exerted by a given ambient flow on a motionless and a freely 
moving particle, respectively. The matrix elements of these 
operators are the single particle scattering coefficients Ai a , 
with I = 1, 2, 3, 4, ... and a = 1, 2, 3 JHf, see Appendix|A] 
for the explicit relation. If the ambient fluid flow is decom- 
posed into the multipole components 12411 . labeled by er, / and 
m = 0, ±1, then Ai a determine the corresponding mul- 
tipoles of the fluid velocity, reflected (scattered) by a parti- 
cle immersed in a given ambient flow. This is why A\ a are 
called "scattering coefficients" Ir9l l23ll . In the multipole ap- 
proach, differences in the internal structure of particles (e.g. 
solid, liquid, gas, porous, core-shell, stick-slip), are fully ac- 
counted by different scattering coefficients. The other parts of 
the multipole algorithm need not to be changed. In particular, 
the scattering coefficeints A10, An and A20 determine the 
single-particle translational, Dq, and rotational, Dq, diffusion 
coefficients and the intrinsic viscosity, [57], by the following 
relations, 



3 M g 



(7) 



D* 
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k B T^ , (8) 
£>5 = fc B T/4 (9) 

with the corresponding single-particle mobility coefficients, 

871-770^20 
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(10) 
(11) 
(12) 
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For core-shell particles, the dependence of all the scattering 
coefficeints Ai a on (x, a, b) have been determined analytically 
inRef. H. 



III. SINGLE CORE-SHELL PARTICLES WITH THIN 
SHELLS 

In this section, we perform asymptotic expansion of the 
core-shell model in the thin-shell limit, i.e. when 



a,K = const and e — > 0. 



(13) 
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The asymptotic expansion of all the single-particle scatter- 
ing coefficients Ai cr (x,a,b), derived in Ref. [9], in the thin- 
shell limit ( fT3l . is the following, 



A la (x,a,b) = A^(a)[l + a la x 2 e 3 /3 + 
B l2 (x,a,b) = B n s {a)[l + p l2 x 2 e 3 /3 + . 



where 1 = 1,2, a = 0,1,2, 



an 

Pa 



21 - 1, 

a n = 21 + 1, 

21 + 3, 



(14) 
(15) 



(16) 
(17) 
(18) 



and the terms omitted in Eqs. (fl4l-(fT5Tl scale as ~ x 2 e 4 . 

Above, A^(a) and B n s {a) denote the scattering coeffi- 
cients for a hard sphere of radius a i23ll . listed in Appendix 
lAl Taking into account that 



We find out that all the scattering coefficients Ai a (x, a, b), 
up to terms ~ (x 2 e 3 ), are determined by the same effective ra- 
dius d25l l. In particular, also the hydrodynamic radius defined 
by the single-particle rotational diffusion coefficient fllll . 
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1/3 

11 ■ 



(28) 



As discussed in the previous section, the many-body hydro- 
dynamic interactions depend on the particle internal structure 
only through the scattering coeffcients. The result is that in 
the thin-shell expansion ( fl~3l the hydrodynamic interactions 
between core-shell particles consist of the dominant contri- 
bution from the hard spheres with the effective radius, plus a 
small correction which depends on a single parameter x 2 e 4 . 

In the next section, the single-particle scattering coeffi- 
cients will be used to construct a model for the concentrated 
suspensions of core-shell particles. 





- a 21 - 1 


= a ai °, 


(19) 


Aft(a),A%(a) - 


- a 2l+1 


= a°" 1 =a ai \ 


(20) 




- a 2l+3 


= a ai \ 


(21) 



and using the relation, 

o*(l + *^)*^° [a{l+^)] k , 



(22) 



with tp = x 2 e 3 /3 and k = ai a ,(3i a , we can combine the 
leading terms of the asymptotic expansion (I14M151 I. and write 
it in the following compact form, 



A la (x, a, b) 
Bi 2 (x,a,b) 



A%(a eS )[l + 0(e% 
(a eff )[l + 0(e% 



r>hs 
Bl2 



with the universal a e ff for all the multipole indices, 



deff/ a - 



l + lx 2 e 3 



0{e% 



(23) 
(24) 



(25) 



and the C(e 4 ) terms ~ x 2 e 4 . The above relations show that 
in the leading order of the thin-shell expansion, the core-shell 
scattering coefficients are equal to the scattering coefficients 
Abides) and _B^ s (a e ff) of a hard sphere with the effective 
radius a e s- This result is the important generalization of the 
thin-shell analysis performed in Ref. [8] for the intrinsic vis- 
cosity and translational self-diffusion. 

In Ref. d8|] it was shown that the hydrodynamic radii, de- 
fined by the intrinsic viscosity, 
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5 



1/3 



(26) 



and the translational diffusion coefficient of single core-shell 
particles ifToll . 



a eff,t — g^iOj 



(27) 



have the same leading terms d25T l in the thin-shell expansion. 



IV. ANNULUS MODEL 

For semi-dilute and concentrated suspensions, direct inter- 
actions between particles have to be taken into account. In the 
annulus model B25I1 . these interactions are taken into account 
by introducing a non-overlap radius b, which surrounds a hard 
sphere of radius c, which is smaller than b. The shell of the 
relative thickness 5, 



(29) 



is filled with the same fluid as outside the annulus particle, 
see Fig. [2] The annulus model can in general be applied for 
particles of a different internal structure and a different type 
of direct interactions. 




FIG. 2: The annulus particle: the solid core of radius c is surrounded 
by a completely permeable (ft = 0) shell with the no-overlap radius b. 

Transport coefficients of semi-dilute and concentrated core- 
shell suspensions, such as the translational & rotational self- 
diffusion, sedimentation and viscosity coefficients, are evalu- 
ated as statistical averages of the corresponding many-particle 
mobilities. The leading contributions to the moblilities of the 
core-shell particles with the outer and inner radii b and a (see 
Fig. Q]) come from the hydrodynamic infractions of the hard- 
spheres with effective radii a e s, larger than a and smaller than 
b. On the other hand, the averaging takes into account the no- 
overlap condition - the particle centers have to be separated 
by a distance at least twice as large as the outer shell radius 
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b. Therefore, the dynamics of many core-shell particles is ap- 
proximated within the (effective) annulus model with c = a e s 
and the non-overlap radius b the same as for the core-shell 
particle, see Fig. [3] 



The explicit theoretical expressions for the two-particle 
virial coefficients A = X t , Xk , X r and A,, are derived in Ap- 
pendix [B] In Sec. IV CI they are used to evaluate X(x, e) nu- 
merically in a wide range of the parameters. 





core-shell: a, b, K 



annulus: a^f, b, K=0 



FIG. 3: In the dynamics of concentrated suspensions, a core-shell 
particle with a thin shell is well-approximated by the annulus model 
with the hydrodynamic inner radius c = a e ft. 

The corresponding effective hydrodynamic shell thickness, 

b - a eff 



e e ff : 



a e ff 



(30) 



V. EXAMPLE OF THE THIN-SHELL ANALYSIS : VIRAL 
EXPANSION OF THE TRANSPORT COEFFICIENTS 

In this Section, we will explicitly perform the thin-shell ex- 
pansion for the two-particle virial contribution to the trans- 
lational and rotatinal self-diffusion, sedimentation and high- 
frequency viscosity, and compare the results with the effective 
annulus approximation. 



A. Virial expansion 

The virial expansion of the transport coefficients for sus- 
pensions of core-shell particles is carried out using the inner 
radius a as the reference one, i.e. we expand in powers of the 
volume fraction, 



-nna 



(31) 



The short-time translational and rotational self-diffusion 
coefficients, D t and D r , the sedimentation coefficient K, and 
the high-frequency viscosity ijoo have the following virial ex- 
pansion, 



Dt/Dl 
K 

Voo/V0 



1 + X t( f> a + 0(4>1), (32) 

l + X r <t> a + 0{4> 2 a ), (33) 

l + X K <Pa + 0(4> 2 a ), (34) 

1 + Wa + \<t>l + 0(4>l), (35) 



with the single core-shell particle translational & rotational 
diffusion coefficients, and the intrinsic viscosity defined in 
Eqs. d8l>- ([T2l i. The bar over the virial coefficients reminds that 
they correspond to the expansion with respect to volume frac- 
tion based on the inner radius (rather than the outer one). 



B. Annulus model 

An annulus particle, introduced in Sec. [IV] and illustrated 
in Fig |2 is described by two parameters: the hard core ra- 
dius c and the no-overlap radius b. The corresponding shell 
thickness <5 is defined in Eq. j29t . With the use of the volume 
fraction <fi c = Anne 3 / '3, the annulus two-particle virial coeffi- 
cients X A (S) of the translational and rotational self-diffusion, 
sedimentation and viscosity are defined by the analogs of 
Eqs. (l32l-(l3"5T> with <fi a — > <fi c , and the single core-shell par- 
ticle coefficients (l8l>-([T2li replaced by the corresponding val- 
ues for a hard sphere with radius c. The two-particle virial 
coefficients X A (S) are evaluated and listed in table U 



TABLE I: The annulus model: Two-particle contributions to the 
short-time translational and rotational self-diffusion, X A and X A , 
sedimentation, X^, and high-frequency viscosity, A^. 



Xi 



xi 



A r 



X'n 



0.000 
0.020 
0.040 
0.060 
0.080 
0.100 
0.120 
0.140 
0.160 
0.180 
0.200 
0.220 
0.240 
0.260 
0.280 
0.300 
0.320 
0.340 
0.360 
0.380 
0.400 
0.420 
0.440 
0.460 
0.480 
0.500 



1.8315 
1.7821 
1.7397 
1.7015 
1.6664 
1.6337 
1.6030 
1.5740 
1.5465 
1.5203 
1.4952 
1.4712 
1.4481 
1.4259 
1.4045 
1.3838 
1.3638 
1.3445 
1.3257 
1.3075 
1.2898 
1.2727 
1.2560 
1.2397 
1.2239 
1.2085 



-6.5464 

-6.7672 

-6.9943 

-7.2271 

-7.4654 

-7.7090 

-7.9578 

-8.2118 

-8.4710 

-8.7353 

-9.0048 

-9.2793 

-9.5589 

-9.8436 

•10.1334 

•10.4282 

■10.7281 

•11.0331 

•11.3430 

■11.6580 

■11.9780 

■12.3031 

■12.6332 

■12.9683 

■13.3084 

■13.6535 



■0.6305 
■0.5483 
0.4916 
■0.4460 
■0.4079 
■0.3752 
■0.3466 
■0.3214 
■0.2990 
■0.2789 
■0.2608 
■0.2444 
■0.2294 
■0.2158 
■0.2033 
■0.1919 
■0.1813 
■0.1716 
■0.1626 
■0.1542 
■0.1465 
■0.1392 
■0.1325 
■0.1262 
■0.1203 
■0.1148 



5.0021 
4.8123 
4.6554 
4.5187 
4.3972 
4.2877 
4.1881 
4.0970 
4.0132 
3.9357 
3.8639 
3.7970 
3.7347 
3.6765 
3.6220 
3.5708 
3.5228 
3.4777 
3.4352 
3.3952 
3.3574 
3.3217 
3.2880 
3.2560 
3.2258 
3.1972 



The effective annulus intrinsic viscosity is the same as for 
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FIG. 4: The two-particle virial coefficients X(x,e), as functions of the shell thickness e, for suspensions of core-shell particles with x = 
(solid line), and x = 3, 10, 20, 50, 100 (symbols), in comparison to the effective annulus model, A^ f (dashed lines). Top: the translational and 
rotational self-diffusion coefficients, \± and A r . Bottom: the sedimentation and high-frequency viscosity coefficients, Xk and A,,. 



a solid sphere, 

[V\ A (S) = |. (36) 

As discussed in Sec. UlTland illustrated in Fig. [3] the effec- 
tive annulus approximation corresponds to c = a e ft, with the 
effective radius given by Eqs. d25ll-(l28T>. 

To compare the core-shell virial coefficients X(x, e) with 
the corresponding values for the effective annulus model, we 
first evaluate the effective radius a e ff and the corresponding ef- 
fective shell thickness e e ft from Eqs. (I26l>-(l28ll and d30l >. then 
we evaluate from Table U the two-particle virial coefficients 
A (eeff) for the annulus model, and finally perform the ap- 
propriate rescaling of the volume fraction QcIf = ima^ a /3 
to (j) a = Anna 3 /3 (which is used in defining X(x, e) for the 
core-shell model), 

A e 1 f ee A^(e eff )(^) a , (37) 



with a — 3 for At, Xk, A r and a = 6 for A^. The effective 
annulus intrinsic viscosity, 




C. The core-shell particles with thin shells 

Based on the previous sections, we now conclude that for 
core-shell particles with thin shells, the virial coefficients [rj\ 
and A are well-approximated by the corresponding values fol- 
lowing from the effective annulus model, 

[vMv]ts and A«A& fore^O, (39) 

with the neglected terms scaling as ~ (x 2 e 4 ). With the use 
of the expressions from Appendix [B] the two-particle virial 
coefficients X(x, e) have been evaluated. In Figs.|4]and|5] the 
values of A for the core-shell particles are compared with their 
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FIG. 5: Unlargement of Fig.[4]for very thin shells. 



analogs A^ f for the effective annulus model. It is clear that 
the effective annulus approximation is very accurate, with the 
range of validity even wider than it could be expected from 
the thin-shell expansion ( fT3l l. 



of two-particle interactions. Large differences between the 
accurate effective annulus approximation to A r , and the cor- 
responding predictions of the imprecise effective hard-sphere 
model are illustrated in Fig. [6] 



VI. CONCLUSIONS 

There are two main conlcusions following from this work. 
The first one is that for core-shell particles with thin shells, 
characterized by large permeabilities, the effective annulus 
model very accurately approximates the self-diffusion, sedi- 
mentation and viscosity coefficients of core-shell suspensions. 
This general statement has been derived from the asymptotic 
expansion of the scattering coefficients in the thin-shell limit 
(a, k — const, e — > 0), and explicitly illustrated by evalu- 
ation of the two-particle virial contributions to the transport 
coefficients, X(x, e), which are shown in Fig. [4] Actually, 
the range of validity of the effective annulus approximation, 
A = A^ f [l + 0(e 4 )], with the neglected terms ~ x 2 e 4 , ex- 
tends beyond the case of x 2 e 4 < 1. 

The second important outcome is that the effective annu- 
lus model for suspensions of core-shell particles with highly 
permeable shells is much more accurate than just the effec- 
tive hard-core approximation, in which the particle shells are 
totally neglected (as completely permeable and overlapping 
with each other). This effect is explicitly shown on the level 



x=3 

»*" ~" - - - 




- , x=20 

N \x=20 v \ 


— effective annulus model 
— effective hard-sphere model 





0.05 0.1 0.15 0.2 0.25 0.3 



£ 

FIG. 6: For very permeable shells, the effective hard-sphere model 
is not accurate, in contrast to the effective annulus model, compare 
with the right upper panel of Fig. [4] 
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TABLE II: Relative accuracy AA of a virial coefficient A in the effec- 
tive annulus and effective hard sphere models, for e=0.3 and x=10. 















AA t 


AA A - 


AA r 


AA ?) 


effective annulus model 


1% 


0.2% 


3% 


2% 


effective hard sphere model 


14% 


13% 


66% 


19% 



The comparison of the accuracy of both models applied to 
X(x—10, e=0.3) is performed in Table ITU For x < 10 and e > 
0.3, the effective hard-core model is not adequate, especially 
for the rotational self-diffusion. 

The main result of this paper is to derive, in the thin-shell 
limit (a, n = const, e — > 0), validity of the effective annulus 
approximation to the core-shell systems, and to demonstrate 
how to use table Q] of the annulus values to determine accu- 
rately the core-shell transport coefficients in the semi-dilute 



regime (and for relatively thin shells), knowing the geometri- 
cal and the hydrodynamical radii of a core shell-particle. 
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Appendix A: Scattering coefficients 

The single-particle friction operator Zo is expressed in 
terms of the scattering coefficients by the following relation, 



Z (ilma,il'm'a') = 6wS mrn >r)(2l + l) 



21(21-1) 
l + l 



A 



10 







(2l-l)(2l+l)A l2 



(2l-l)(2l+l)A l2 



l(l+l)A n 

(m)(2m) 2 (2z+3) 



o 



21 



Bl2 



r 



(Al) 



and 



Zq — Zq — Zo /x Zo, 



(A2) 



with the only non-zero matrix elements of fi for I = 1 and 

(7 = 0,1, 

Ho(ilm<j, ilm'a) — \Z^(i\raa, ilm'a)] . (A3) 

The explicit (lengthy) expressions for the scattering coeffi- 
cients Ai J7 (x, a, b) for core-shell particles have been derived 
in Ref. [9]. Their asymptotic expansion in the thin-shell limit 
(fT3T i is given in Sec. [HI] Below, we list the scattering coeffi- 
cients for hard spheres of radius a, 



A - 2l + 1 n 2l-l 
Aiq — — - — a 



An = 



2 

2l + 'i 



An — a , 



-a 



21 + 1 



Bl2 — 



21 + 1 



.21+3 



(A4) 



Appendix B: Two-particle virial coefficients 

In this Appendix, we derive expressions for the two-particle 
virial coefficients X(x, e) of the core-shell suspensions. In the 
limit of k=0 (i.e. x=0), the annulus model with the inner ra- 
dius c = a is recovered. 



We introduce dimensionless interparticle distance, 
R=r/2a. The pair distribution function has the form, 



9o(R) 



for R < 1 + e, 

1 for R > 1 + e. 



(Bl) 



which corresponds to the no-overlap condition at the larger 
shell radius b = a(l + e). 

The two-particle virial coefficients are evaluated from the 
following expressions, 



r+oo 

A t (>,e) = 8 / J t (R)R 2 dR, 

Jl+e 

2 4 o 

Xk(x, e) = ^A 12 - -Ai (l + eY 

f + 00 



r+cc 

8 / J K (R)R 2 dR, 
Ji+t 



+ 

r+oo 

X r (x,e) = 8 / J r (R)R 2 dR, 

Jl+e 

2 - - r°° 

X v (x,e) = -[ryf + 24[ry] / J v (R)R 2 dR, 

O Jl+e 



(B2) 

(B3) 
(B4) 
(B5) 
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where 



Jk(R) 



1 



J t (i?) = ^ F Tr M t 1 f ) (R), 
Mo 



Mo 
1 



j txif ) (R)+ / i!r(R)-T (R 



tt(2). 



J r (i2) = ^Tr^(R), 
Mo 



(B6) 
, (B7) 
(B8) 
(B9) 



with the single particle dipole mobility /1q given by Eqs. (0 
and ( [Tol l. In the above expressions, Tr denotes the trace oper- 
ation, and 



T (R) 



1+RR 

8irr]oR ' 



(BIO) 



is the Oseen tensor, with R = R/R. All the mobility coef- 
ficients, the associated functions J and the scattering coeffi- 
cients Ai a in Eqs. ( !B21 i-( lB9t refer to the core-shell particles. 
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